Influence of electronic correlations on orbital polarizations in the parent and doped iron pnictides 
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Orbital polarization and electronic correlation are two essential aspects in understanding the normal state and 
superconducting properties of multiorbital FeAs-based superconductors. In this Letter, we present a systematical 
study on the orbital polarization of iron pnictides from weak to strong Coulomb correlations within the KRSB 
approach. The magnetic phase diagram of the two-orbital model for LaFeAsO clearly shows that the striped 
antiferromagnetic metallic phase with considerable orbital polarization exists over a wide doping range. The 
reversal of the orbital polarization occurs in the intermediate correlation regime.A small crystal field splitting 
enhances the orbital density wave order. 



o 

o 
i 

s : 

a ■ 



c 

o 
o 



> 

oo 
m 

o 



X 
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Since the discovery of high-r c superconductivity (SC) in 
doped iron pnictides 13 , the iron-based high-T f SC material 
has been the subject of intensive researches in the past two 
years. Among various normal-state and SC properties of iron 
pnictides, the roles of the electronic correlation and the orbital 
polarization have been the central topics. The parent phases 
of cuprates are Heisenberg antiferromagnets (AFM) and insu- 
lators, and the important role of the electronic correlation and 
dominant single-orbital hole with the x 2 - y 2 symmetry were 
well established and verified. Unlike the high- T c cuprates, the 
parent phase of iron pnictide SC are a bad metal 2 , and usu- 
ally exhibits spin-density-wave (SDW)-type or striped AFM 
ordering 3,5 . What is the role of the electronic correlation and 
how many Fe 3d orbits are involved in the iron pnictides are 
two essential and key aspects to understand various properties 
in the SDW metallic and SC phase. 

Up to date, experimental evidences on the role of electronic 
correlation in parent phase and SC phases are controversial. 
Indeed, the metallic conductivity^, small magnetic moment 
of Fe spins 5 , X-ray absorption spectroscopy (XAS) and reso- 
nant inelastic X-ray scattering (RIXS) of Fe 4 in undoped iron 
pnictides seem to suggest a weak correlation between Fe 3d 
electrons. From the resonant X-ray emission spectroscopy 
(RXES) of Fe 3d Z-2,3 edges, Kurmaev et al. 6 suggested that 
the 111 phase is a weakly or at most moderately correlated 
system. However, many more experimental clues demon- 
strated the band narrowing effect, suggesting the important 
roles of the 3d electronic correlation 4,6,7 . In fact, in addi- 
tion to the band narrowing effect, the existence of interme- 
diate magnetic moment in the 1111 , 122 and 111 phases also 
shows the importance of the electronic correlation. There- 
fore one should employ an approach applicable for both weak 
and strong Coulomb interactions so as to uncover the roles of 
electronic correlations on the electronic states and magnetic 
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configuration in iron pnictides. 

On the other hand, it is still under hot debate in literature 
on how many orbits are involved and what the orbital symme- 
try of the bands near Ep is in the parent and doped iron pnic- 
tides. Theoretically, two-orbital 8 , three-orbital 9 , four-orbital 10 
and even five-orbital 11,12 tight-binding model were proposed. 
Most of the multi-orbital tight-binding models captured the 
major Fermi surface and dominant band structures characters 
near £> of the 1111 and 122 phases. To distinguish the valid- 
ity of these models is essential to reveal the symmetry char- 
acter and the number of the orbits involved in various iron 
pnictides. 

The role of the electronic correlation in the groundstate 
properties of parent iron pnictides was discussed in very 
recent paper a 14 ' 15 . To uncover the roles of the electronic 
correlation and orbital polarization in high-Tc iron pnictide 
SC, we apply the Kotliar-Ruckenstein's slave boson (KRSB) 
method 13 on the two-orbital Hubbard models. The KRSB ap- 
proach and its extension may provide a useful tool to treat the 
multi-orbital Hubbard model over a wide correlation range. 
It also has a few of advantages in dealing with multiple or- 
bitals, nontrivial magnetic configurations and spatial correla- 
tions. Our recent numerical ansatz developed for the multi- 
orbital KRSB solution overcomes the convergency of many 
parameters in minimizing the groundstate energy, suitable for 
treating arbitrary coulomb interaction in the presence of in- 
terorbital hoppings or hybridizations and various magnetic 
configurations. 

In this paper, we first present the magnetic phase diagram 
of two-orbital model and demonstrate that in the intermediate 
electronic correlation, the striped AFM (SAFM) phase with 
small magnetic moment is stable against the paramagnetic 
(PM), ferromagnetic (FM) and Neel AFM phases, accompa- 
nied with a small orbital polarization, or a ferro-orbital density 
wave (F-ODW) order. A small crystal field splitting arising 
from the lattice distortion further stabilizes the SAFM and F- 
ODW ground state. The evolution of the orbital occupation 
and Fermi surface with the electron correlation shows that the 
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contribution of the d xz orbital symmetry is larger than that of 
the d yz orbital. These results clearly show the importance of 
the intermediate electronic correlation in the parent and doped 
phases of iron pnictides. 

We start with two-orbital model Hamiltonian H — Ho + Hi 
applied for iron pnictides in real space with 



<i J>apo~ 



Hi = U ^ j Wtof/tfal, + V (V - JH^mT^niao-riifrj' 
ia iao-'{a>(i) 
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where Cj creates an electron with orbital index a and spin 
cr at lattice site i, n laa - is the corresponding occupation num- 
ber operator. The hopping integral for the orbital a and (3 is 
denoted by t a p. The intraband (inter-band) Coulomb repul- 
sion and Hund's rule coupling are denoted by U (£/') and Jh , 
respectively. Here we set U' — U — Uh- 

To reflect the multi-orbital character of iron pnictides, we 
extend the single-orbital KRSB approach 13 to the two-orbital 
Hubbard models for various magnetic configurations. In the 
multiorbital Hubbard model, a few of auxiliary boson field 
operators representing the possibilities of various electron oc- 
cupations are introduced, such as e, p, d, b, t, q, which denote 
the possibilities of none, single, double, triplicate, quaternity 
occupations. With these auxiliary boson fields, an original 
fermion operator can be expressed as: 



Q 



(1 - QiatrTzfL 



(3) 



Here fJ is the new slaved fermion operator and Q iaiT is an 
auxiliary particle number operator 13 . Projecting the original 
fermion operators into these boson field and fermion field op- 
erators, one could not only obtain an effective Hamiltonian, 
but also get the groundstate energy in the saddle point ap- 
proximation with the normalization condition and the fermion 
number constraints 13 . Here we employ a generalized La- 
grange multiplier method to enforce these constraint condi- 
tions, thus the interorbital hoppings and crystal field splitting 
can be treated on the same foot. The fermion occupation num- 
ber is constrained with the penalty function method. To en- 
force the normalization condition, we have a boundary con- 
strained condition: 

acr a acrcr' acr 

We use optimizing method to get the minimized groundstate 
energy. To search for possible ground state, we study four 
different magnetic configurations, including the PM or non- 
magnetic, FM, Neel AFM and SAFM phases so as to find the 
stable ground state. 
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FIG. 1: Magnetic phase diagram as a function of particle density n 
and correlated U for J H = 0.2U. PM, SAFM(LM) and SAFM(HM) 
stand for paramagnetic metallic, striped antiferromagnetic metallic 
phases with magnetic moment m < l/u B and m > l/u B , respectively. 



As the penalty factor is sufficiently large, the particle num- 
ber constraint is satisfied. Nevertheless, the problem of 
boundary constrained condition is still difficult. In numeri- 
cally searching the global minima of the groundstate energy, 
we employ the pattern search method, together with the gra- 
dient method and the Rosenbrock method. If the optimiz- 
ing point is on the boundary, we move one step inward the 
high-dimensional ellipsoid and the equipotential plane. Since 
the first axis of the new local orthogonal coordinate system 
of Rosenbrock method directs to the negative gradient direc- 
tion, the corresponding algorithm is simple. To investigate the 
groundstate electronic and magnetic properties of iron pnic- 
tides, we adopt the hopping parameters of Zhang et al^ for 
the two-orbital situation in the present slave boson scheme. In 
the two-orbital situation, the orbits 1 and 2 refer to the d x z and 
d y z components, respectively. 

We first obtain the magnetic-phase diagram of the two- 
orbital system through comparing the energies of the PM (or 
nonmagnetic), FM, Neel AFM and SAFM configurations. In 
the present symmetric orbital situation, the penalty factors 
adopted for the two degenerate orbits are identical. This will 
save a great deal of computation time in minimizing the total 
energy. 

FigtUdisplays the zero-temperature phase diagram of U vs 
n for Zhang et a/.'s 8 hopping parameters. We consider the 
full Hund's coupling Jh, which equals to Q.2U . One notices 
that only the PM metallic phase is stable when the on-site 
Coulomb correlation U < U c \ « 1.2 eV. When U > U c \, 
the SAFM phases are stable over wide doping and interaction 
ranges. This SAFM region can be divided into high magnetic 
moment one (SAFM(HM)) with m > lyUg and low magnetic 
moment one (SAFM(LM)) with m < 1/-Ib- The latter lying 
in U c i < U < U C 2 is the most interesting and marked with a 
shadow since the magnetic moments of the parent and doped 
phases of the most of iron pnictides, including the 1111, 122 
and 111 phases, fall into this region. In the present phase di- 
agram, it is also shown that the SAFM(LM) regime is limited 
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FIG. 2: Electronic correlation dependence of total energies of PM 
phase (dashed line) and SAFM phase (solid line) and that of the en- 
ergy difference {E pm -Esafm) between these two magnetic configu- 
rations at n = 1.9, 1.95,2,2.05 and 2.10, respectively. 



to a very narrow range. 

Especial interest is found around the particle filling number 
of n = 2. One notices that the phase diagram is not symmet- 
rical about the half filling (n = 2) due to the presence of the 
next nearest-neighbor hopping. When U > U c \, the system 
enters the spin gapped insulating phase. In this insulating re- 
gion the gap between the spin-down and spin-up subbands, the 
sublattice magnetic moment, the spin exchange splitting and 
the quasiparticle bandwidth increase with the increase of U. 
When the electronic correlation becomes so strong that U > 5 
eV, the system enters a full gapped region, very similar to the 
conventional Mott insulating phase. In the insulating region, 
with the further increase of the Coulomb interaction U, the 
spin exchange splitting approaches a constant and sublattice 
magnetic moment is saturated. In comparison with the con- 
ventional PM Mott phase, the quasiparticle bandwidth in this 
SAFM Mott phase gradually decrease with the increase of U 
again. As shown in the following, the orbital polarization is 
observed in a wide doping and an interaction regime, while 
no orbital polarization is seen at n = 2 when the correlation is 
strong. 

For more clarity, we plot the correlation dependence of the 
groundstate energies in the PM structure and the SAFM one 
in Fig|2] Also an energy difference, E PM - E SAFM , between 
these two phases is displayed. One finds that when the elec- 
tronic correlation is weak, or U < U c j, the energy differ- 
ence is negative and very small. So the PM phase is stable 
against the SAFM one. When the correlation reaches to the 
intermediate strength that U > U c i , the energy difference be- 
comes positive, the SAFM phase becomes more stable than 
the PM phase. Furthermore, one could estimate the magnitude 
of the magnetic couplings between iron spins. The magnetic 
coupling strength crucially depends on the doping concentra- 
tion and electron correlations, as seen in Figj2] For a typical 
parent compound LaOFeAs, the on-site Coulomb correlation 
between Fe 3d electrons is estimated to 2 ~ 3 eV— . This 
gives rise to the magnetic coupling strength about 20 meV 
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FIG. 3: Sublattice magnetic moment in the SAFM phase as a func- 
tion of the Coulomb repulsion strength U. Theoretical parameters 
are the same to Fig.l. 



at n = 2, in the same magnitude order of the spin coupling 
between Fe ions in LaOFeAs deduced from the neutron scat- 
tering experiment 5 . 

Fig[3] shows the evolution of sublattice magnetic moments 
with the increase of the electronic correlation in the SAFM 
phase. It shows that the system remains PM without any mag- 
netic ordering unless U becomes larger than U c u which in- 
creases with hole doping. In the realistic parameter range, i.e., 
intermediate correlation regime, the magnetization decreases 
with the increase of hole doping, in agreement with the ex- 
periment observation. We notice that the minimal value of 
U c i is about 1.2 eV at n = 2, which is considerably larger 
than the mean-field approximation results for the two-orbital 
model 1 ^^. It is also found that the sublattice magnetic mo- 
ment increases continuously until U reaches U C 2 = 2.4 eV, 
where it changes discontinuously. This behavior is also found 
by Yu et al. in the mean-field solution in the four-orbital 
model 17 . Further, the upper critical value for the SAFM 
regime with m < l/xg is 2.4 eV for n = 2.0, significantly 
larger than that obtained by the mean-field approximation^ 7 *^, 
suggesting that the quantum fluctuations are well considered 
within the present KRSB approach. At n — 2.0, the sublattice 
magnetic moment is larger than 1 fig when U > U C 2 = 2.4 
eV. The SAFM system enters the high magnetic moment re- 
gion, this is also shown in FigQ] By monitoring the density 
of states (DOS) near the Fermi surface, we find that due to 
strong correlation, the system becomes insulating at n — 2.0. 
We notice such a fact that when the systems cross over from 
the SAFM insulating phase to the full-gap Mott one, the sub- 
lattice magnetic moments of the systems with n — 1.9 and 
2.1, or n = 1.95 and 2.05, approach the same values, as seen 
in Fig. 3, showing that the electron-hole symmetry restores in 
the strong correlation regime. 

In the present two-orbital model, we can address the or- 
bital polarization or orbital ordering problem in the parent 
and doped phases of LaOFeAs. We find that in the weakly 
correlated PM phase , due to the equivalence of the present 
two orbitals, the symmetry of two orbitals is not broken, no 
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FIG. 4: Magnetic and crystal field dependence of orbital polariza- 
tions or itinerant orbital density wave ordering in two-orbital models. 
Theoretical parameters are the same to Fig.l. 
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FIG. 5: The DOS at the Fermi surface as a function of U at different 
doping concentrations. 



orbital polarization or obrital ordering is observed. As soon 
as the system enters the intermediate correlated SAFM phase, 
the orbital polarization is observed, as seen in Fig.4. Fig. 4(a) 
shows that in the SAFM phase, the orbital polarization consid- 
erably depends on the Coulomb correlation and particle fill- 
ing. When n < 2, the orbital polarization is negative, or the 
d yz component is enhanced with the increase of U. At « = 2, 
the orbital polarization only occurs in the intermediate corre- 



lation regime, and the polarization is reversed at f/ Q = 1 -9 eV; 
when Ud < U < Uq the polarization is dominated by the d xz 
component, while it is dominated by the d yz component when 
Uo < U < U C 2- When U > U C 2, due to strong correlation, the 
two electrons occupy two orbitals with the same possibility, 
giving rise to a vanishing polarization. This behavior was also 
found in the conventional mean-field approximation 19 . When 
n > 2, we notice that the polarization also reverses the sign at 
about Uq. Though such an orbital polarization of ordering is 
weak, it is greatly enhanced in the presence of a small crystal 
field splitting, as shown in Fig. 4(b). From Fig.4(b) one finds 
that the orbital polarization do not reverse over a wide corre- 
lation regime in the presence of crystal field splitting, while a 
kink at U = still exists for various splitting. 

These results clearly demonstrate that the orbital polariza- 
tion or ferro-orbital ordering is accompanied with the onset of 
the SAFM ordering in the intermediate correlated iron pnic- 
tides. Since in the intermediate correlation regime, the system 
is metallic, a question naturally arises that what the essential 
of the F-ODW ordering is in metallic iron pnictides. Obvi- 
ously, such an orbital polarization is essentially itinerant, as 
the itinerant FM in transition metals. Hence the orbital or- 
dering in iron pnictides is characteristic of density wave. As 
we pointed out in an early paper 20 , such an itinerant ferro- 
orbital ordering can occur as soon as a Stoner-like condition 
is satisfied 20 . Besides, in Fig|5]we also display the correlation 
dependence of the quasiparticle DOS near E F in the present 
two-orbital systems. It is found that the DOS increases slowly 
in the weak correlated PM region. With the increase of U, the 
DOS first steeply decreases and then slightly increases in the 
metallic SAFM ordered phase with small magnetic moment, 
demonstrating a pseudogap behavior. The pseudogap in the 
intermediate correlation regime turns to a full gap in the spin- 
density-wave gap insulating phase when U > 2.7 eV. From 
the analysis on the spin-dependent partial DOS, we find this 
full gap is triggered by spin polarization. When U > 4.1 eV, 
the systen enters a conventional Mott insulator phase. 

In summary, we have shown that in the parent and doped 
phases of iron pnictides, the intermediate electronic cor- 
relation favors a striped antiferromagnetic configuration 
accompanied with a small structural distortion and weak 
ferro-orbital ordering. The orbital polarization is the charac- 
teristics of a ferro-orbital density wave. 
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